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CN , Abstract. I present some new results, from simulations and observa- 

?H ' tions, of the constraints bar pattern speeds place on the dark matter 

content of disk galaxies. 
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Introduction 



Near infra-red images reveal bars in over 50% of disk galaxies (Knapen et al. 
2000; Eskridge et al. 2000). The principal dynamical quantity for barred (SB) 
galaxies is the pattern speed of the bar, Q p , usually parametrized by the ratio 
TZ = Di/as (where Dl is the corotation radius and as is the semi-major axis 
| of the bar). A bar is termed fast when 1.0 <1Z^ 1.4, while, for a larger value of 

TZ, a bar is said to be slow. All measurements of bar pattern speeds have found 
fast bars (e.g. Gerssen 2002), a result which has been interpreted as evidence 
for maximum disks in SB galaxies (Debattista & Sellwood 1998, 2000). 



2. Bar Dynamical Friction in iV-Body Simulations 



Fig. 1 plots the evolution of Dl and a# for a near-maximum disk, and demon- 
strates that fast bars can survive for a Hubble time in such systems. This 
simulation used a hybrid grid code, allowing us to achieve much higher resolu- 
tion than in previous work. At the end of the simulation, which lasted some 44 
initial rotation periods of the bar, 1Z = 1.3 ± 0.1. Although the bar slows down 
by ~ 35% in this case, secondary bar growth can still preserve a fast bar in this 
case. But secondary bar growth has its limitations: for a flat rotation curve, 
Dl\/Dl2 = ^p 2 /^p 15 so a strongly braked bar would need to grow ever longer, 
which, however, makes friction even stronger (Debattista & Sellwood 2000). 

Weinberg &; Katz (2001) argued that TV-body simulations require N p ^> 
4 x 10 6 particles to resolve correctly the bar-halo interaction, while simulations 
at lower N p are largely noise-driven. In Fig. 2 I plot the fractional drop in Q p 
between early and late times in a series of simulations in which only N p was 
varied. No significant change in the evolution of Q p is seen, despite the factor 
of 100 difference in N p , suggesting that Weinberg & Katz (2001) have been 
too pessimistic in their assessment of noise effects. This can be understood by 
recognizing that a bar slows continuously, so that resonances are broadened. If 
resonances had been sharp, then indeed large A p 's would be required, but for 
broad resonances, more modest JV p 's suffice. 
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Figure 1. Evolution of as (open circles with error bars) and Dl 
(filled circles) for a massive disk simulation. At t = 1000, 1Z = 1.3±0.1. 
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Figure 2. The change in U p , due to halo dynamical friction, between 
early and late times for a series of simulations in which N p has been 
varied by a factor of 100. 



3. A Post- Interact ion Galaxy 



Following Ostriker h Peebles (1973), it is sometimes thought that unbarred (SA) 
galaxies must be stabilized by massive dark matter (DM) halos; however, bars 
can also be suppressed by rapidly rising rotation curves (Toomre 1981; Sellwood 
& Evans 2001). Noguchi (1987) showed that disks stabilized by massive DM 
halos can still form bars via interactions, and that such bars would be slow. Thus 
SB galaxies with evidence of past interactions provide a means to statistically 
test whether SA galaxies are stabilized by massive halos. 

NGC 1023 shows evidence of a past interaction (Sancisi et al. 1984), without 
being at present significantly perturbed. Debattista et al. (2002) measured 
lip for this galaxy, using the method of Tremaine & Weinberg (1984). The 
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Tremaine- Weinberg (TW) method is contained in the following simple equation: 

. . _ fZoHY)Vio S (X,Y)Z(X,Y) dXdY ^ V 
pSm '" C 00 h(Y)XE(X,Y)dX dY ~X [ ) 

where V\ os is the line-of-sight velocity, £ is the surface brightness, h{Y) is an 
arbitrary weight function, i is the inclination and (X, Y) are galaxy-centered 
coordinates along the disk's major and minor axes respectively. We obtained 
3 slit spectra parallel to the disk's major axis, for each of which we computed 
V and X. Then, plotting V versus X gives a straight line with slope Op sin i = 
4.7 ± 1.7 km s _1 arcsec _1 . From the rotation curve, corrected for asymmetric 
drift, this gives Dl = 531^° . Multi-band photometry gives ob = 69" ± 5"; 

therefore 1Z = 0.77t°;25, consistent with a fast bar. Thus NGC 1023 must have 
a maximum disk; moreover, if the bar formed in the interaction (which cannot 
be ascertained), then it cannot have been stabilized by a massive DM halo. 



4. Pattern Speed in the Milky Way 

As shown by Kuijken & Tremaine (1994), the TW method can also be applied to 
the Milky Way Galaxy (MWG). For discrete tracers, the TW method becomes: 

AV = n p R -v LSR = ^ 

J2if{n)Vr,i 
J2if(ri) sin k cos bi 

where Rq is the Sun-MWG center distance, Vlsr is the tangential velocity of 
the local standard of rest (LSR), ulsr is the radial velocity of the LSR, f(n) 
is the observational detection probability (which need not be known), v T) i is 
the heliocentric radial velocity of a discrete tracer, and (k,bi) are its Galactic 
coordinates. A survey which satisfies the condition of uniform sampling, as 
required by the method, is the ATCA/VLA OH 1612 MHz survey (Sevenster et 
al. 1997a,b & 2001), covering |/| < 45° and \b\ < 3°. The properties of OH/IR 
stars make them well suited for Galactic studies. They have bright, maser 
emission at 1612.23 MHz, which is insensitive to interstellar extinction; the 
double-peaked profile of the emission permits easy identification; and they are 
concentrated in the inner MWG, as a result of the Galactic metallicity gradient. 

We extracted from the ATCA/VLA OH 1612 MHz survey a sample of ~ 250 
OH/IR stars which are older than 0.8 Gyr and have a flux density higher than 
0.16 Jy; these criteria give OH/IR stars between ~ 4 and ~ 10 kpc away from 
the Sun. Applying Eqn. 2 to this sample, as shown in Fig. 3, we obtained 
Ay = 252 ± 41 km s _1 , where the error was estimated by resampling. If we 
then assume Vlsr/-Ro = 220/8 km s _1 kpc -1 (from SgrA* motion, Backer et al. 
1999; Reid et al. 1999 and Cepheid proper motions, Feast & Whitelock 1997) 
and mlsr = (from SgrA* HI absorption spectrum, Radhakrishnan et al. 1980), 
we obtain Q p = 59 ± 5 ± 10 km s _1 kpc -1 , where the last error is our estimated 
systematic error. This £l p is consistent with the hydro dynamical simulations of 
Bissantz et al. (2002), which also require 72. = 1.0 ± 0.1. 
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Figure 3. The TW analysis of the OH/IR stars for changing Iq (the 
maximum \l\ in Eqn. 2). From left to right are V, /C, the number of 
OH/IR stars and the resulting KjV 
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